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ABSTRACT 

This paper describes the position sensorless 
algorithms presently used in the motor control for the 
NASA "in-house" development work of the flywheel energy 
storage system. At zero and low speeds a signal 
injection technique, the self-sensing method, is used to 
determine rotor position. At higher speeds, an open loop 
estimate of the back EMF of the machine is made to 
determine the rotor position. At start up, the rotor is set 
to a known position by commanding dc into one of the 
phase windings. Experimental results up to 52,000 rpm 
are presented. 

INTRODUCTION 

NASA is in the process of developing the technology 
necessary to make flywheels attractive for use as energy 
storage devices in space applications. One of the 
necessary components is the motor/generator and its 
associated control (Kenny, 2001). The permanent 
magnet (PM) machine is an attractive choice for the 
motor/generator because of its high power density and 
efficiency. The PM machine can be accurately controlled 
with a fast response and low torque ripple through the 
use of the vector control algorithm (field orientation) 
(Krishnan,1999). However, field orientation requires rotor 
position information with a resolution on the order of a few 
degrees (Jahns, 1994). Position feedback with this 
resolution for PM machines is typically provided by using 
either an encoder or a resolver. However there are speed 
limitations on these devices and even in speed ranges 
where position feedback devices are feasible, a 
sensorless approach still has the advantage of reducing 
the complexity of the system by eliminating a hardware 
component and its associated failure modes, volume and 
mass. 


A variety of sensorless techniques have been 
proposed in the literature and it is still an active field of 
research (Johnson, 1999). This work combines two of 
the techniques to achieve sensorless control at all 
speeds: the self-sensing technique for low and zero 
speed operation (Jansen, 1995) and a back EMF based 
method for higher speeds. A similar combination is 
described in (Patel, 2001). 

NOMENCLATURE 

L q is the q-axis machine inductance, henries. 

L<j is the d-axis machine inductance, henries. 
f* y represents a q-axis scalar quantity where f can be 
voltage, v, current, i, or flux, X; x can be the rotor 
reference frame, r, or the stator reference frame, s, or 
the negative sequence reference frame, neg; and y is 
either the stator variable, s, or the rotor variable, r. 

f| y represents a d-axis scalar quantity where f can be 
voltage, v, current, i, or flux, X; x can be the rotor 
reference frame, r, or the stator reference frame, s, or 
the negative sequence reference frame, neg; and y is 
either the stator variable, s, or the rotor variable, r. 

fq dy represents the d-q vector quantity where f can be 
voltage, v, current, i, or flux, X; x can be the rotor 
reference frame, r, or the stator reference frame, s, or 
the negative sequence reference frame, neg; and y is 
either the stator variable, s, or the rotor variable, r. 

A superscript asterisk, *, represents a commanded 
value. 

A carrot above a symbol, A , represents an estimated 
quantity. 
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Vqds c is the commanded high frequency carrier signal 
voltage vector for the self-sensing technique, 
p is the derivative operator, d/dt. 

0 r is the angle between the stator q-axis and the rotor 
q- axis, radians. 

X a f is the flux linkage due to the rotor magnets, volt-sec 
% is the electrical rotor speed, radians/second. 
co c is the carrier signal frequency, radians/second. 
co e is the fundamental electrical excitation frequency, 
radians/second. 


is independent of the fundamental operation of the 
machine and is thus present at all speeds. The carrier 
signal, V q( j s c , is added in a feedforward fashion to the 
commanded fundamental voltage and both are 
synthesized by the PWM action of the inverter switches. 
The block diagram of the overall system is shown in Fig. 
1. The vector quantities shown in Fig. 1 are defined by 
(1) where f can represent current, flux or voltage. Both 
scalar and vector forms are used in this paper. 

fqds = fqs " jfds ("I) 


SELF-SENSING TECHNIQUE 

The self-sensing method of position estimation in AC 
machines was first described in (Jansen, 1995). It is a 
sensorless technique that relies on a magnetic saliency 
present in the machine to provide position-dependent 
information which can be extracted with the appropriate 
signal processing. Unlike the back EMF techniques, the 
method does not rely on fundamental excitation of the 
machine so position estimation is possible at low and 
zero speeds. Flowever, the method does require an 
additional signal into the machine which will lead to 
additional losses. 

A necessary requirement of the self-sensing 
technique is a known magnetic saliency in the machine 
so that a position (or flux) dependent inductance results. 
Most permanent magnet machines inherently have some 
magnetic saliency with interior permanent magnet 
machines exhibiting the most saliency (L q = 2 to 3 times 
Ld) and surface mount permanent magnet machines 
exhibiting a very small saliency (L q =1.05 to 1.15 times L d ) 
(Krishnan,1999). The machine used in this work is a 4 
pole surface mount permanent magnet machine with Lq = 
139 pH and L d = 116 pH. The larger the saliency, the 
easier it is to track the position-dependent signal. 
However, even with the small saliency in this machine it 
is possible to successfully estimate the rotor position. 

The other necessary requirement of the self-sensing 
technique is a high frequency carrier signal in addition to 
the fundamental excitation. The carrier excitation signal 


A high frequency current with a position-dependent 
phase variation results from the injection of the carrier 
signal voltage and the position-dependent inductance 
variation in the machine. The expression for this high 
frequency current will now be derived. 

The machine inductance vector, expressed in the 
rotor reference frame, is given in (2). 



The expression for the stator flux linkage at the carrier 
frequency is shown in (3). 




Lq 0 iq S 

OlJL ids- 


(3) 


The transformation between the rotor reference frame and 
the stator reference frame is shown in (4) 
(Krishnan,1999). 


f q r sl = [ cos(< 

- f d r J L-sin( 


I cos(0 r ) sin(0 r ) f qs 

.- sin(0 r ) cos(0 r )J Lf ds J 

Substituting (4) into (3) transforms the equations from the 
rotor reference frame into stator frame as follows: 
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Figure 1: SYSTEM BLOCK 
DIAGRAM 






and 


cos(0 r ) sin(0 r ) 

^qsl f Lq 0 If cos(0 r ) s 

m(0 r )fiq S l 

. -sin(0 r ) cos(0 r )_ 

>4] L 0 L d J_-sin(0 r ) cos(0 r )j[j.s j 

kslj cos(0 r ) -s 

in(0 r ) L q oTcos(0 r )s 

in(0r)Tiqsl 

UjL si "( 0 r) c < 

)s(0 r ) _ _ 0 L d JL-sin(0 r ) 

OS(0r)JLisJ 


Thus the stationary frame expression for the inductance 
is given in (7) and shows the position dependent 
variation. 


Vm £L qd 

° >C ( £L q d) 2 - ( ALq,}) 2 


| V m AL qd 

° >C ( 2L qd) 2 - ( AL qd) 2 

Equation (11) can be expressed in vector form by using 
the definition given in (1) as shown in (12). 

iqds c = -jlcp^ - jlcnej (20r ■® ct) (12) 


s J cos(0 r ) -sin(0 r )TL q 0 T cos(0 r ) sin(0 r )l 
c l d X sin(0 r ) cos(0 r ) JL 0 L,jJ_-sin(0 r ) cos(0 r )J 

SL qd + AL qd cos( 20 r ) -AL qd sin( 20 r ) 

-AL qd sin( 20 r ) ZL qd - AL qd cos( 20 r ) 

Ln + L d Ln - L d 

where IL qd =^-, AL qd =^- (7) 

The machine model in the stator reference frame at the 
carrier frequency is given in (8). 

Vqds c = R s iqds_c + L qd piqds_c (8) 

Substituting (7) into (8) and neglecting the stator 
resistance (the stator resistance is small relative to the 
inductive impedance at the carrier frequency) gives the 
relationship between the high frequency injected carrier 
voltage and the resulting high frequency current. 


Vqs_c 

ZL qd + AL qd cos( 20 r ) -AL qd sin( 20 r ) 

-AL qd sin( 20 r ) SL qd - AL qd cos( 20 r ) 

The applied high frequency carrier signal is a balanced 
polyphase voltage given by (10) where V m is the peak 
value of the applied signal. 



V qs _c 


V m cos( m c t) 

-Vdsc- 


-V m sin( co c t) 


Substituting (10) into (9) and solving for the high 
frequency current results in the following. 


'qs_c 

l cp sin(m c t) + l cn sin(-co c t + 20 r ) 

-'ds_c-l 

l cp cos(co c t) + l cn cos( -ra c t + 20 r ) 


It can be seen from (11) and (12) that there are two 
components to the high frequency current. Both 
components are at the carrier signal frequency, eo c . The 
l cp term is the "positive sequence" component, that is, it 
is rotating in the same direction as the applied voltage 
(+<» c ). The l cn term is the "negative sequence" 
component, that is, it is rotating in the opposite direction 
of the applied voltage (-<o c ). The position information, 0 r , 
is seen to be contained in the negative sequence 
component. 

Equation (12) represents the high frequency current 
in the machine due to the application of the carrier signal 
voltage and the presence of the magnetic saliency. The 
total stator current will also contain a fundamental 
component with an excitation frequency co e as shown in 
(13). 

iqds = - jlcpe^ 1 - jlcnej (20r ' ^ 3 ) 

The filtering of this current will be described next. 

The positive sequence component of the current, 
•jIcp eJC ° ct » contains no position information nor is it 
necessary for the fundamental operation of the machine. 
A synchronous frame filter is used to effectively remove 
this component as shown in Fig. 2 ("HPF" is a high pass 
filter). 

I 6 e j| CD eW - jl cn e + K 2e r ~ m c t) 
l e e + K“et-®c‘) _ j, j| cn efj(20 r -2to c t) 


| e e +j(fo e t-oa c t) . j| cne +j(20 r -2to c t) 


Figure 2: Positive sequence synchronous frame filter 

The remaining current is seen to consist of the 
fundamental component and the negative sequence 

component: l e ei® et - jlcn 6 ^ 201 " Wct ^ The fundamental 
component is necessary to regulate the fundamental 
operation of the machine. The negative sequence 
component is necessary to estimate the rotor position. 
Thus two additional filtering steps are needed as shown in 
Figure 3. To extract the fundamental current, a 
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Figure 3: Negative sequence synchronous frame filter 

AND FUNDAMENTAL SYNCHRONOUS FRAME FILTER 

negative sequence synchronous frame filter is used. To 
extract the negative sequence current, a fundamental 
synchronous frame filter is used. 

The position information is extracted from the 
negative sequence current by using a heterodyning 
process to create an error signal which is then used in an 
observer structure as shown in Figure 4. In a flywheel 
system, where the majority of the electromechanical 
torque, T em , is used to accelerate or decelerate the 
flywheel, the use of the estimated torque in the observer 
improves the position and speed estimate by feeding 
forward the acceleration information to eliminate the 
position and speed estimate lag. The observer structure 
shown in Figure 4 can also be implemented without the 
torque feedforward term which results in essentially a 
phase locked loop structure. Note that a speed estimate 
is also available in this structure and is used as feedback 
in the outer speed control loop of the motor. 



SPEED OBSERVER 

The scalar representation of the cross product shown in 
Fig. 4 is given in (14) where (0 r - 0 r ) is assumed to be 
small. 

^-(-jl cn ej (20r) X - jl cn eK 2 0r>) = 

.neg ^ .neg 

-j7*-cos(20 r ) - ■4—sin(20 r ) = sin 2(0 r - 0f) ~ 2(0 r - 0 r ) (14) 


START-UP PROCEDURE 

The rotor position angle, 0 r , is defined as the angle 
between the stationary frame q-axis and the rotor frame 
q-axis. The self-sensing method determines 0 r based on 
the inductance difference between the d- and q-axes. The 
method does not inherently differentiate between the 
positive q-axis and the negative q-axis. This can be 
seen from (14) where the phase of the filtered current is 
proportional to twice the rotor position. This means that a 
rotor position angle of 0 electrical degrees will result in 
the same negative sequence current as a rotor position 
angle of 180°. If the rotor position is at 0 electrical 
degrees and the self-sensing estimate locks onto 180°, 
the result will be that the machine will spin in the negative 
direction for a positive speed command. So an initiation 
test must be done at start up to determine if the self¬ 
sensing estimate has locked onto the correct angle. 

The initiation test is done as follows. First the rotor 
is moved to a known position. This is done by 
commanding a small dc voltage onto the stator windings 
which effectively sets up a dc magnetic field in the air 
gap. The rotor magnets will align with this field and thus 
the rotor will move to a known position (Matsui, 1996). 
By commanding a small additional duty cycle to the 
positive rail switch on phase A of the machine, a dc field 
is established in the direction of the stator +q-axis. The 
rotor magnets align with this field which results in a 90° 
angle between the stator q-axis and the rotor q-axis (ie. 
0 r = 90°). The self-sensing algorithm will estimate either 
0 r = 90° or 0 r = -90° for the rotor in this position. If the 
estimate is equal to -90°, then 180° is added to the self¬ 
sensing estimate and that result is used in the control 
algorithm for 0 r . If the estimate is equal to 90°, then 
there is no offset angle added to the estimate before 
being used in the control algorithm. 

BACK EMF TECHNIQUE 

The back EMF technique is based on the fact that 
the spinning rotor magnets will generate a voltage (the 
back EMF) at the terminals of the machine. This voltage 
can be integrated to find the stator flux (Wu, 1991). From 
the stator flux vector and an estimate of the torque 
angle, the rotor flux position can be estimated (Patel, 
2001). This will be described next. 

The stator reference frame equation describing the 
relationship between the stator voltage, current and flux 
linkages for the fundamental excitation of the machine is 
given in (15). 

v qds = iqds^s + P^-qds (15) 

The stator flux can be estimated by integrating the stator 
voltage less the IR drop as shown in (16). 

^qds v qds ' iqds^s^t (16) 

In the simplest implementation (used in this work), 
the stator voltage, v q< i s , is assumed to be equal to v q< i s , 
the commanded voltage, which is known in the controller. 
This is equivalent to assuming a constant dc bus voltage 
and no voltage loss in the inverter. Improvements in the 
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flux estimate by using either better estimates of stator 
voltage or by actually measuring stator voltage are 
possible and described in the literature. 

Additionally, the open loop integration described by 
(16) was actually implemented as a low pass filter with a 
50 Hz bandwidth. This was found experimentally to result 
in a reliable stator flux estimate at speeds above about 
500 rpm. 

However, it is knowledge of the position of the rotor 
flux, X a f, which is necessary for proper field orientation. 
This can be found by expressing the stator flux in the 
rotor reference frame wherein the d-axis component of 
the stator flux vector is aligned with the rotor flux, X a f. 



Figure 6: Speed observer 

also used as an input into a speed observer which is 
similar in structure to Fig. 4 and shown in Fig. 6. Again, a 
torque estimate can be used to enable zero lag properties 
of the speed estimate (Lorenz, 1991). 


4 = L q'qs ( 17 ) 

^-ds = Ultls ' ^-af (18) 

The angle between the stator flux vector, Xqds, and X a f is 
the torque angle which is found from (17) and (18). 



These relationships can be best shown graphically as 
in Fig. 5. For field orientation control, the d-axis of the 
rotor reference frame is defined to be coincident with the 
rotor flux. To transform from the stator frame to the rotor 
reference frame, the rotor angle, 0 r , must be known. This 
is the angle between the stator frame d-axis and the rotor 
frame d-axis (or, equivalently, the angle between the 
stator frame q-axis and the rotor frame q-axis) as shown 
in Fig. 5. Mathematically, this angle can be found from 
(20) as shown. 


0 r 



(20) 


The estimated rotor angle, 0 r , is used for the 
reference frame transformation as shown in Fig. 1 and is 



Figure 5: Vector diagram showing flux, current, angles 

AND REFERENCE FRAME AXES 


EXPERIMENTAL RESULTS 

The sensorless algorithms were implemented to 
control this motor precisely because position feedback 
was not available. Thus it was not possible to 
independently verify the exact accuracy of the position 
estimates. However, a general indication of the fidelity of 
the estimates is presented here by observing the motor 
operation and the position estimates under certain 
conditions, including successful operation to 52,000 rpm. 

One method to gauge accuracy is to observe the 
position estimates while the machine is freely 
decelerating ("spinning down"). Under this condition, the 
fundamental stator current is commanded to zero and 
(15) reduces to: 




r x af cos(0 r ) 

) r X af sin(0 r ). 


(21) 


Thus there is a way to relate the position estimate to 
the fundamental phase voltage applied to the machine. 
The actual phase voltages were not measured. However 
the commanded voltage, Vq S , is available in the controller 
and was assumed to be approximately equal to the actual 
phase voltage, V an . 

Figure 7 shows the rotor position estimates during a 
spin down test from 1500 rpm. It can be seen that there 
is a phase difference between the two estimates. A 

400 1 - 1 - 1 - T - T - 1 



Figure 7: 0 r estimates for spin down from 1500 rpm. 
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lime, seconds 

Figure 8: Expected and actual q-axis phase voltages 
DURING SPIN DOWN FROM 1500 RPM 

comparison of the phase shift between the actual rotor 
position and the estimates can be found by using (21) 
and the estimated rotor angles to calculate the predicted 
Vq, The predicted voltages and V^, s are plotted in Fig. 8. 
It can be seen that the back EMF method estimate leads 
the actual by about 44° while the self-sensing lags it by 
about 24°. 

The phase lead in the back EMF method is due to the 
implementation of (16). A low pass filter with a bandwidth 
of about 50 Hz is used instead of the pure integration 
described in (16). The dc offsets in the measurements 
cause erroneous results if a pure integration is used. A 
pure integrator has a phase shift of 90° however the low 
pass filter approaches a 90° phase shift only at 
frequencies above the bandwidth. Thus the back EMF 
estimate improves with higher speeds. This can be seen 
in Figure 9 where the data from a spin down from 7500 
rpm is shown and the phase lead is reduced to 
approximately 10°. 

The self-sensing error can be attributed to two main 
factors. First, the resistance of the machine was 
neglected in the derivation of (11) and (12) and this will 
lead to a slight phase lag in the estimate. Second, the 
filtering of the current to determine the position estimate 



Figure 9: Expected and actual q-axis phase voltages 
DURING SPIN DOWN FROM 7500 RPM 



Figure 10: Phase current Figure 11: Phase current 
DURING TRANSITION WITH NO DURING TRANSITION WITH 


OFFSET IN 9,- OFFSET IN 0 r 

introduces phase errors. The relationship between these 
effects and the position estimate error is deterministic 
and the estimate could be improved with additional tuning 
and testing. However, in this application exact accuracy 
was not important because the self-sensing method was 
used only at low speeds. 

The position and speed feedback are changed from 
the self-sensing estimate to the back EMF estimate 
manually, typically at 1500 rpm. Fig. 10 shows the phase 
current during the transition with no adjustment for the 
phase errors of the estimates. Fig. 11 shows the 
transition with both estimates adjusted appropriately. It 
can be seen that there is a spike in the current without 
the adjustments but the effect is small and it is not 
necessary to adjust for the phase errors to successfully 
transition between the methods. 

The speed estimate for both techniques is shown in 
Fig. 12 at a constant speed of 1500 rpm. Both speed 
estimates are within +/- 1% of the commanded value. 
The accuracy of the speed estimate can be best seen by 
examining a frequency spectrum of the current. In a 
synchronous 4 pole machine, the fundamental frequency 
occurs at twice the mechanical frequency. Thus at 1500 
rpm, the expected fundamental current frequency is 50 
Hz. Fig. 13 shows this to be the case. The second 
harmonic oscillations in the back EMF speed estimate 
decrease significantly at higher speeds, probably due to 
the better position estimate. The bandwidth of the speed 
loop is relatively low so the speed controller does not 
respond to the oscillations in the speed estimate. 



time, seconds 

Figure 12: Estimated rotor speed for a commanded 
VALUE OF 1 500 RPM 
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Figure 13: Phase current Figure 14: Phase current 
FREQUENCY SPECTRUM, LOW FREQUENCY SPECTRUM, HIGH 
FREQUENCY RANGE FREQUENCY RANGE 


The amount of current necessary for the self-sensing 
technique can be seen from the high frequency end of 
the spectrum as shown in Fig. 14 and is approximately 
3.2 amps. This results in a negative sequence current of 
approximately .25 amps which is the input to the observer 
as shown in Fig. 4. After the feedback is switched to the 
back EMF estimates, the high frequency voltage is turned 
off and the only fundamental current remains. 

Results from operation at 52,000 rpm are presented 
in Figs. 15 through 17. The fundamental frequency is 
approximately 1733 Hz. These results were achieved 
with a 290 volt dc bus, 40 kHz controller sample rate and 
65 kHz switching frequency using a MOSFET based 
inverter. 



0 0.002 0.004 0.006 0.008 0.01 

time, seconds 

Figure 15: Phase current for 52,000 rpm operation 


I 

I 

I 



Figure 16: Estimated 
ROTOR SPEED FOR A 
COMMANDED VALUE OF 
52,000 RPM 


!' 5 


Figure 17: Phase current 
FREQUENCY SPECTRUM FOR 
OPERATION AT 52,000 RPM 



CONCLUSIONS 

NASA Glenn Research Center has demonstrated 
hybrid sensorless control of a surface-mount PM machine 
for a flywheel application at speeds up to 52,000 rpm in 
an in-house test facility. The PM machine is started and 
accelerated to 1,500 rpm using the self-sensing 
technique of position estimation in a closed loop speed 
control with field orientation. Then the feedback is 
changed to the back EMF estimates and the additional 
high frequency voltage necessary for the self-sensing 
method is turned off. The machine is then controlled to 
the desired speed using the back EMF estimates. The 
self-sensing method is initially calibrated by commanding 
the rotor to a known position and determining if an offset 
angle of 180° is necessary. This hybrid technique allows 
the machine to be operated over its full speed range, 
from start up to full speed, without a sensor for position 
feedback. This was critical in this application as this 
motor was not provided with a resolver or an encoder. 

Future work will include automating the transfer from 
self-sensing feedback to back EMF feedback at a pre¬ 
determined speed. In addition, the exact phase errors for 
the different methods will be studied more thoroughly 
once a machine with a resolver is available. 
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